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Abstract 

Inspired by recent experiments by Geim et al. we discuss the classical the- 
ory of the Hall effect of a 2 dimensional electron gas in an inhomogeneous 
magnetic field. The field modulation is in the form of flux tubes created by 
a superconductor overlayer. We find that an approach, where the vortices 
are treated as individual scatterers contributing to the collision term in the 
Boltzman equation will not work — it leads to a vanishing Hall constant at 
T = 0. If the field is treated as a smooth contribution to the driving term 
in the Boltzman equation, the classical Hall constant emerges, in agreement 
with experiments when the Fermi wavelength is short in comparison with all 

other lengths in the problem. 
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I. INTRODUCTION 



We will in this paper discuss the Hall effect in the case, where the applied magnetic 
field is spatially varying. The theory will be based on the Boltzman equation, and we shall 
only consider the 2D case. Our motivation comes from recent experiments by Andrei Geim 
et al.0, who put a superconductor over a 2 dimensional electron gas, and then measured 
the Hall voltage. The superconductor will only allow the magnetic field to penetrate in 
Abrikosov vortices, thereby modulating the field. For high fields the vortices are strongly 
overlapping and the field is only slowly varying, so that the usual Hall coefficient is to be 
expected. This is indeed what is seen experimentally and it is the result of our calculations. 
For low fields (below 100 G), where the vortices start to become spatially separated, the Hall 
coefficient depends on the 2D electron density and therefore on the de Broglie wavelength 
of the electrons at the Fermi surface. When the de Broglie wavelength is comparable to 
or greater than the diameter of a vortex, the Hall effect is reduced by an almost field 
independent fraction in low fields. The fraction is about 80% in the electron gas with the 
smallest density experimentally obtainable. Since the effect depends on the electron density 
one might expect that it is a quantum mechanical effect. The only way to put quantum 
mechanics into the Boltzmann equation is through the scattering cross sections. Khaetskiii 
has proposed treating the spatially separated vortices as asymmetric scatterers. At de 
Broglie wavelengths much shorter than the vortex diameter the electron will be scattered 
asymmetrically and in accordance with a classical picture. At the opposite limit, first treated 
by Aharonov and Bohmi, where the wavelength is much larger than a vortex diameter the 
scattering is symmetric. Khaetskii and earlier Kuptsov and Moiseevll showed that the degree 
of asymmetry gradually disappears as the diameter of the vortex is reduced in comparison 
with the electron wavelength. Khaetskii's idea is that this reduced asymmetric scattering 
can account for the reduced Hall effect. His calculations show that this is indeed possible in 
a classical gas. We have extended his calculations to a degenerate gas obeying Fermi-Dirac 
statistics. Here the Pauli principle, and its restrictions on the scattering, will reduce the 
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calculated Hall coefficient to about k-^T/eF, in strong disagreement with the experiment. 

We have subsequently considered the case where the magnetic field is a slowly modulated 
field. This amounts to treating the 5-field as a driving force on the left hand side of the 
Boltzmann equation. Here we find the classical Hall effect corresponding to a homogeneous 
field at all field strenghts and electron densities. This is to be expected since we recover the 
experimental results when the electron wavelength is shorter than a vortex and we henceforth 
can talk about the magnetic field at the electron's position with some confidence, while the 
procedure fails when the electron wavelength is longer than the field modulations. 

II. VORTICES AS SCATTERING CENTERS 

First we will describe the magnetic flux tubes as independent scattering centers, much 
like usual impurities. This of course is supposed to apply only at the low field limit, where 
the tubes are sufficiently far apart. The characteristic feature of scattering off flux tubes is 
that the scattering probability is asymmetric: There is an enhanced probability of electrons 
being scattered to the left. We will model this by an asymmetric scattering probability 
w(k,ip), where tp is the scattering angle and k is the length of the momentum vector. The 
Boltzman equation, linearized in the external electric field, has the familiar form 



= (df\ 



where f° is the equibrilium distribution function. The collision term consists of two parts, a 
flux tube part and a usual impurity scattering part, which we will treat in the relaxation time 
approximation. Denoting the distribution function as f(k,<ft), <p being the angle between k 
and the external field E, we have 
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Here p is the density of fluxtubes, i.e. p = (BA)/((p /2)/ A = eB/h. The important 
difference to the work by Khaetskiii, is the inclusion of the Pauli principle. 



We will solve the equation by Fourier transforming in the angle 0. Introducing 

/»(*) = r t- <!>)> w ^ = r t- e<B * w (*' ®> ^ 

Jo In JO 2,11 

the Boltzman equation (P has the form 
df° 

-evE-L-{8 n>1 + Vi)/ 2 = -ifn(k) - <W )A (4) 
- p((wo(k) - w n (k))f n (k) +J2{w m (k) -W- m (k))f n - m (k)f m (k)j. 
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The functions w n (k) satisfy w n (k)* = W- n (k). Accordingly there are both real and imag- 
inary contributions to the effective relaxation time, due to the flux tubes. The imaginary 
parts have, as pointed out by Khaetskii, a simple interpretation, namely as an effective 
homogeneous magnetic field. Indeed in the Fourier transformed Boltzman equation with a 
homogeneous magnetic field, the magnetic field term has the form: 

e(v x — -imu c f n (k). (5) 

op 

Upon linearization of the last term in @ we get fo = f° and of the other terms only f\ and 
/_i are non- vanishing and they become 

fl{k) = f - l(kr = 21+ip(2f (k)-l)hz( Wl (k))T(e) E > (6) 

where 

r(e)- 1 = p(l-Re( Wl (A;))) + r- 1 . (7) 



It is now straightforward to work out the conductivities. We get 
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and 
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where the bracket (■) is defined by 



and a(e) = pr(e)Im(w 1 (e)). Khaetskiii has shown that in the classical limit a(e F ) = ^ C T - 
The important difference between the result in equation (^) and Khaetskii's result is the 
factor 2/o (e) — 1, which is zero at the Fermi level. This means that the Hall voltage will 
disappear at T = 0. The factor comes from the proper implementation of the Pauli principle. 
In the low-T (T <C e F ) limit, where we approximate a and r by their value at the Fermi 
level, we simply get 

ne 2 r(e F ) Arctan(a(e F )) 

o xx = r~\ — • C 11 ) 

m a(e Fj ) 

By neglecting the term ((2/ (e) — l)a(e)) 2 in the denominator in (§) we get 

k B T ne 2 r(e F ) . . , . 

cr xy < — —a e F . 12 

e F m 

In the experiment by Andrei Geim et al.0 the temperature was 1.3 K. If we use an effective 
mass of 0.07m e , k^I '/e F is less than 0.1. The electron mobilities were in the range of 40-100 
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and the magnetic field was swept from G to 200 G. Consequently a is in the range of 0- 
2. In Figure 1 we have plotted the Hall resistance (normalized to the classical value B e g/ne) 
as a function of the magnetic field. From here it is seen that the Hall effect is reduced by 
a factor of about ksT/ep. To illustrate the crossover to the non-degenerate electron gas 
case treated by Khaetskii we have plotted the same quantity as a function of temperature 
in Figure 2. 

The conclusion is that asymmetric scattering does not give rise to a Hall effect in a 
degenerate electron gas. This result is not in agreement with experiments. To explain this 
result we take as a simple model the scattering to the left through the same angle </>o at each 
scattering event. 

w{$) oc 5(0 - O ) (13) 

The effect of the — eE field is to make more electrons go in it's direction (6 = 0). The effect 
of the scattering and therefore of the magnetic field is to oppose this effect by scattering 
electrons out of the 8 = direction. In particular (considering only magnetic scattering) 



'Of 
dt 



< (14) 

mag-coll 



From (|2|) we have with our model scattering (|T^) that 

(§0 (*, 0) = f(k, -0 O )(1 - /(*, 0)) - f(k, 0)(1 - f(k, O )) (15) 

\ / coll-mag 

We want to determine the angle dependence of f(k,9). If f{k,9) is small {k > k Fermi ) it is 
the angle dependence of the f(k, 9) outside the parentheses that dominates and we therefore 
drop the parentheses and get from ( |14|) 



f{k,-<f> )-f{k,0)<0 (16) 
f(k,-(j> )<f(k,0). (17) 

Consequently the electrons have a tendency to move to the left. This classical picture is 
due to the fact that the parentheses we have neglected are exactly the contribution from the 
Pauli principle. If, on the other hand, f(k, 9) is close to 1 (k < kpermi) the Pauli contribution 
dominates and we consequently drop the prefactor to the parentheses. 

(i-/(M))-(i-/(Mo))<o (is) 

f(k,<f> o )<f(k,0) (19) 

The electrons now tend to move to the right. The reason is that in order to scatter in a 
dense Fermi gas it is essential that there are few electrons a scattering angle away. The 
electrons above and below the Fermi surface thus move in opposite directions (the 1 — 2fo(e) 
factor) and the net asymmetry is very small (it arises solely from the difference in speed 
below and above the Fermi surface) — the Hall effect has disappeared. 

III. MAGNETIC FIELD AS A DRIVING FORCE 

From the above section we conclude that it is not correct to treat the vortices as scattering 
centers, that discontinously changes a electron's position in phasespace. In this section we 
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are going to discuss a complementary approach, where we treat the inhomogeneous field as 
a driving force, that changes a electron's position continously in phase space. The approach 
is totally classical and is supposed only to apply in a dense electron gas, where the electron 
de Broglie wavelength is short in comparison with the length over which the magnetic field 
varies. We will assume that the magnetic field is random, correlated over lengths comparable 
to the effective London length. 

To the usual linear order in the electric field and in the deviation g(f, k) from equilibrium 
f°(f, k), the Boltzmann equation is 

f |(f, I) - • ^(fc) - e(t7 x B(r)) • |(f, fc) ^ ||j (f,k). (20) 

The collision contribution we will treat as scattering against fixed impurities. Accordingly 
in polar coordinates in the fc-space 



'df\ 



., , (?,e) = p r^w(cf>)(mo -<!>)- me)), (21) 
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where p is the density of scatterers. We have suppressed the k dependence. 

We will write the magnetic field as B(r) = B° + 5B(f), where B° is the average magnetic 
field. The Boltzmann equation can now be written 

v cos9-^-(f, k) + v sin 9-^-(f, k) + u c -^-(f, k) 
ox ay 69 

eSBdg r df° k flpn fdf\ ^ 

H ^\ r i k ) - ev^— cos 9 E x (r) - ev— — sin 9E y (r) = — (r,k) 

m o9 oe oe y at J „ 

(22) 

with uo c = eB°/m. 

In (^) we introduce the hermitian operator 

D =4 +ir <( c ° se l +aiae -Q- (23) 

where r c = v/u> c is the classical cyclotron radius in the magnetic field B°. We now want 
to simulate the effect of the collisions by a relaxation time approximation, so we add and 
subtract a relaxation time contribution and arrive at 



D + 



9 = X 



i df2 /r . . _. . .5B dg 

—ev^(E x cos9 + E y sm9)-i—^ 



u; r r w r V at 



(24) 



(25) 
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Here r can be chosen arbitrarily, but later the usual value of r will emerge as a natural 
choice. The eigenfunctions of D with eigenvalue n (n integer) are 

1 



exp (iq ■ f — m6> — ir c (fc sin + ir c q y cos 0) , 



(26) 



where A is the area of the electron gas. Therefore g can be written 

g = J d? J d9'G(r,r*,9,9') X (f*,9') 

with the Green function 



(27) 



q,n 
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(28) 



-in(0-9') 
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i -in(e-e') 



where 



Po{4>) = r + r c 



V 



sin 9 
- cos 9 



+ r c 



( - sin (9-<f))^ 



+ COS (0 - 0) 



(29) 



is the classical cyclotron orbit in the homogeneous field B° parametrised by the momentum 
coordinates (angles). We assume uj c t is positive and get with = 9 — 9', using Poisson's 
summation formula 



1 ^ e" 



^ iexp(-J^) 



2ttV?- + — ~exp(-^)-l 
where |0| is the value of in [0, 2tt[. We, finally, have 

r2w 



(30) 



exp (— 2tt/uj c t) — 1 Jo 



d(f>exp (-<I>/u c t)x(P${4>), 9 -<t>) 



(31) 



The physical interpretation of this formula is that you assume that the electrons move 
along their classical trajectory in a homogeneous magnetic field -B . The correction to the 
local electron density is obtained by adding field corrections from the neighbouring points 
according to the number of electrons arriving from neighbouring points to your fieldpoint. 
The prefactor arises because we only integrate around the classical circular orbit once — we 
could drop it and instead integrate to infinity. 

For our later choice of relaxation time the mean free path is very long (at least 2//m and 
normally more than 10/im) so in fact we make a field average along the classical trajectory. 
Now the cyclotron radius r c is of the order 2/im and therefore much bigger than the magnetic 
correlation length £, which is of the order 0.1/im, so the system is strongly selfaveraging. 

We now average and get 



<g> = 



7 Z*^ 71 " -* 

— / d4>e~^ < X {Pe{<t>),0-4>)> 

— 1 JO 



g— 2-k/u} c t 



/•27T $ ( QfQ eV 

- / d<pe «ct i-± (< E x > cos (9 - 0)+ < E v > sin (0 - (f>)) 

1 Jo \ oe uj r 



g — 2tt/u) c t 



with 



+ £l±(e - <t>) + ^ r ^W)(< g{e-4>-4>')>-< g(o - <p) >) 

U! C T u c Jo Air 



# u (f, #,r) = rev-^-cos6 



de \ 1 + (cj c t) 2 1 + (cu c r 
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Notice that g° is the distribution function in a homogeneous magnetic field in the relaxation 
time approximation (with relaxation time r). We now decompose the angle part of the 

— * 

/c-space in Fourier components 

g(6) = J29ne- m9 (33) 

n 

w(^) = 5>„e- in * (34) 



Integrating out <f> and <f>' we get 



<g(0)> =g°(f,e,T) + j:- 1 



ie~ my n < 5Bg n > 



--in B° 

UJ C T 

| < g n > e~ ind I p(w n - w ) + _L 



in 



U) r T 



The n = Fourier component is trivial (< g >=< g >) and for n 7^ we have 

1 



<9n> =9n(r) + — 



= 9° n (r) + - 



— in 



1 



in < 5Bg n > p(w n - 1) + \ 

H 21 < 9n > 



B° 



UJ r . 



in 



in < SBSg n > p(w n - 1) + \ 



B° 



+ 



< 9n > 



where we, in the last step, have used < 5B >= 
The current is determined by g\. 

<-°° kdk r 2n 



Jx + IJy 



o 2ir 2 Jo 
00 kdk 

9i- 



d6 g(k, 6){v cos 9 + iv sin 9) 



f°° kdk 

/ V 

Jo it 



If we choose 



- = — = p(l-w 1 )=p d9(l - cos 9)w{9) 

T Ti JO 



we get 



n/ \ * 5B(r) 

<9i> = 9i(n) + < -^T^9i(r) > . 

LO c Tl 

In the appendix we have calculated the leading term in < 5g\ (r) > to be 

e (-27r/o; c Ti) — I J ^ 



5° 5° 



=< #i > O.O61 



/ 10 15 m -2 



n 



(35) 



(36) 



(37) 
(38) 



(39) 



(40) 



(41) 
(42) 

(43) 



We use here that for randomly distributed gaussian vortices each earring half a flux quantum, 



< S£ ff} S) >= 2^25 £ = O.lyum and that r c = ^^-^Jj^nF^ m ^ ne experiment by Geim 
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et al. J 1Qls ™ 2 varies between 0.5 and 1.7. Therefore, the average electron distribution is 
nearly as in the homogeneous case. If we use the approximation ([11]) in (|4(]) and assume 
that T\ can be treated as a constant in the energy integration in fl38|) we get that 



1 + T— TZw^T / ^-^sin0< p y V "n > (44) 



i_ e (- 2 ^ ri )y r B° B° 

^ = i . r #e <-^>cos0 < > . (45) 

The integral in ( pE5| ) is of the same order of magnitude as the prefactor in (pE3|), implying 
the magnetoresistance is a few procent greater than it would have been in a homogeneous 
magnetic field. Because the ^-integration in (pHD is restricted by the correlation function to 
an interval from to the integal is about (j-^j m 400 times smaller than the prefactor 
in (|43D . The approximation (f|l]) is therefore not the dominant contribution to the devia- 
tions in the Hall effect from the homogeneous case. Consequently, the deviation from the 
homogeneous Hall effect is at most a few promille. This is in perfect agreement with the 
fact that the experiment by Geim et al. showed no deviations from the homogeneous result 
in a dense electron gas. 



IV. CONCLUSION 

In the experiment by Geim et al.0 a reduced Hall effect in the Abrikosov vortex modulated 
field is only observed at electron densities below 4 • 10 15 m -2 and in a magnectic field of less 
than 100 Gauss. That is when the external magnetic field varies appreciably within a de 
Broglie wavelength of the electrons at the Fermi surface. In this regime it is expected that 
the Boltzmann equation description breaks down, but outside this regime our treatment of 
the vortices simply as a modulated magnetic field in the Boltzmann equation agrees with the 
experiment. To explain the reduced Hall effect one has to incorporate some kind of quantum 
mechanics. We have shown that it is not feasable to describe the vortices as scatterers and 
hide all the quantum mechanics in the calculation of the scattering cross sections. 
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A full quantum treatment should certainly include multiple coherent scattering by the 
vortices, because single scattering is contained in the present Boltzman calculation. It is 
also clear that the Hall constant is reduced, since in the limit of very thin vortices (or what 
amounts to the same, a very dilute electron gas) the Hall constant will vanish. In this limit 
the time symmetry breaking will vanish, because one can without any change in the physics 
reverse the direction of the field by placing an infinitely thin Dirac vortex carrying one flux 
quantum h/e at each of the external vortices that carries half a flux quantum; the Dirac 
vortices having a field in the opposite direction of the external field. 

We acknowledge discussions with Mads Brandbyge, Erland Brun Hansen, Ayoe Hoff, 
Dung-Hai Lee, Poul Erik Lindelof, Mads Nielsen and Rafael Taboryski. 



APPENDIX A: 

In this appendix we are going to calculate the correlation function 

SB(r) , 1 5B(r) _ „ N ,„ /A1 , 

To do this we assume the higher order correlation functions factorize the second order 
correlation function out and we henceforth have the gausssian result 

< 5B(r)$(B) >= Jdy< SB(r)SB® >< -^L- > . (A2) 

Using this in (|Al[) we get that 

-£f !^/*<»«»«>><S?>- < A3 > 



>exp {—(p/uj c T) x 



Now we have from (plD that 

8g(r, 9) i r->r. 



SB® exp (— 2tt/lu c t) — l Jo 



SB® C ° S 0) + SB® Sm { ° " 0)) 
. 5{y-Pe{<t>)) dg fi . r ^ n .8B(P B (<j>)) d ( 8g{P B {<j>\0 - 



B° d9 x RVYV ' B° 36 \ 5B® 

- (iJL + 1- (?L 

SB® \ LU C T UJ C \ dt 



colV 
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When the last term is inserted in (|A3|) it is seen as before that if we choose t = T\ the term 
cancels. When we use this expression below we will assume that this kind of cancellation 
can be done, and erase this term. ( |A4|) is an equation to iteratively determine jgnjij- with 
the third term as the driving term. We henceforth expand < - B ^ 5gi (r) > in this term. 
The first order contribution is 



< 



B° exp (— 2ir/u c Ti) — 1 Jo Jo 2n 

<^^>w-A<*>)^<«-*> 

2vr r 2n rjf) , , , 



exp (— 2tt/uj c ti) — 1 Jo Jo 2tt 

5B{r) 5B{P e {4>)) d<g> ro 
<^o B^ > ^^( -^ ^ 

Now < 5B ^ SB ^ ^ > only depends on the distance between r and Pe{4>)- Consequently 
< 6B ^ Q ^ > is independent of 9 and we can move the ^-integral through with the result 
that 

< ^fWi > = -^r\ i r < > ■ <*» 

B" exp (— 2tt/uj c Ti) — 1 Jo B u B" 

< 5B qP sb ^ ^ > is only large within a correlation length £ and it's size is estimated as 

< >• Accordingly, as an order of magnitude estimate we have 

SB(r) s _ £ 5B(f) SB(r) 

< -^T 5 9i(r) >~< 9i>^< ^T^T > ■ ( A7 ) 

This is in our case much less than < gi >. To get the second order contribution we have 
to iterate (Al) once more, putting the driving term back into the first two terms and the 



fourth term on the right hand side of ( |A4| ) . We will first take the fourth term and here we 
get 

< ^Mr) >4<t htcrm = , - ) ; : rd^-^ f dyx 

B" exp (— 2it/uj c Ti) — 1 jo Jo 2tt J 

SB(r) 5B(y) ^ ^ 6B(P e (<j>)) Of 1 /* J ( _^ cn) . 



< 5° B° >< 5° d# l v exp(-27r/cu c r 1 ) - 1 Jo 06 X 
S(y- P 9 {tf> + 0')) |f (Petf + <P'),e-<p- 0')) > (A8) 
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1 ^ > e (-W„) c SB(r)5B(P e ( ( j ) + ( j ) ')) ^ 



(exp (-2tt/o; c ti) - 1) 2 Jo r Jo r 5° 5° 

o 2^ 7i° ^(P e (0 + 0),#-0-0) > . (A9) 

The last integral is, apart from the differentiations, the same as the original integral, just 
now spatially separated. Therefore it is not greater than < g 1 > < SB jP s ^p > - the 
order of magnitude from before. Again the correlation function < 5B (- p e(4'+ < t' )) > j s only 
appreciable within a distance of £. Since both <p and <$' are positive, the two remaining 
integrals are restricted to a region of size Implying that we have the following order of 
magnitude estimate: 

< -jjrfaif) >4th t erm ~ £ [f- < -^r^r >) < * > • (aio) 

If we use that the parenthesis is about 0.1 we get that this term is 10y ~ 10 Q 2 (^ , he. 200 



times smaller than the first order contribution. To take care of the first two terms in (|A4|) 
we use that 

8E x {P e m f d JE x {P 6 {<t>)) 6g (?) (An) 



5B(y) J 6g (z) 5B{y)' 

The last term is treated as above. We find that the first iterate is 0. As explained in the main 
text the first order contribution to < SB ^ 8g\ (r) > mainly influences the magnetoresistance. 
Consequently, higher order terms contribute significiantly to the Hall effect. The most 
important is the first term that arises when you go beyond the gaussian approximation^ : 

< 9i > /** ^ .^(-WOMn)^) < 5B(rl5B(P e (<P))6B(Pe(<P + <P')) ^ 



(exp(^) - l) 2 Jo ^Jo r B° B° B° 

(A12) 



1 In a magnetic field consisting of fluxtubes placed at random 

< 5B(r)5 9l (r) >= £~ x dy n < 5B(r)8B(yi) • • • 5B(y n ) >i-/ to < dB JS% (iJn) >, 

where ./V is the number of fluxes. In this appendix we have calculated the first term in the sum. 
The higher order terms may be calculated in exactly the same manner. 
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The order of magnitude of the relative deviation in the Hall effect from the homogeneous 

2 



case, due to ( |A12| ), is I 0.06 



10 15 m -2 



. Implying that in a dense electron gas the deviation 



in the Hall effect from the homogeneous case is about one promille. 
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FIGURES 

FIG. 1. The theoretical Hall resistivity coming from treating the magnetic fluxtubes as scat- 
terers in the Boltzmann equation, normalized to the classical homogeneous result — as a function 
of a = fiB at a temperature of 0.1, 0.5, 1.0 and 5.0 ep. In the experiments by Andrei Geim et al., 

2 

who used mobilities fx in the range of 40-100 they found in a dense gas the homogeneous result 
at a temperature of less than 0.1 ep. 

FIG. 2. The hall resistivity normalized to the homogeneous value ^ as a function of the 
temperature for a = fiB = 0.5. 
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